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Abstract 



Interpretation of empirical results based on a taxa's lifetime distribution 
shows apparently conflicting results. Species' lifetime is reported to be expo- 
nentially distributed, whereas higher order taxa, such as families or genera, 
follow a broader distribution, compatible with power law decay. We show that 
both these evidences are consistent with a simple evolutionary model that does 
not require specific assumptions on species interaction. The model provides a 
zero-order description of the dynamics of ecological communities and its species 
lifetime distribution can be computed exactly. Different behaviors are found: 
an initial power law, emerging from a random walk type of dynamics, 

which crosses over to a steeper branching process-like regime and finally is 
cutoff by an exponential decay which becomes weaker and weaker as the to- 
tal population increases. Sampling effects can also be taken into account and 
shown to be relevant: if species in the fossil record were sampled according to 
the Fisher log-series distribution, lifetime should be distributed according to a 

power law. Such variability of behaviors in a simple model, combined with 
the scarcity of data available, cast serious doubts on the possibility to validate 
theories of evolution on the basis of species lifetime data. 



Introduction 



Ecosystems have become paradigmatic examples of complex systems, showing organi- 
zation and collective dynamics across very different time and spatial scales [1]. These 
features are captured by non trivial relationships among measurable quantities, which 
take forms familiar to statistical physics. Well known examples include the species- 
area scaling relationship [2,3], allometric relations [4-7], and the occurrence of power 
laws in the distributions of species lifetime and size of extinction events [8-10]. These 
statistical laws have been measured over many orders of magnitude and exhibit sim- 
ilar patterns across very different living ecosystems and also in different quantitative 
studies of fossil records [11]. The ubiquity of these patterns [12] suggests that they 
may be amenable to be studied in a general and a-specific framework. 

In this article, we will address the issue of species (or more general taxa) lifetime 
distribution. Although the analysis of fossil records has recently highlighted several 
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patterns in the evolution of biodiversity, and motivated the proposition of differ- 
ent mechanisms that may have caused these patterns, the functional form of the 
species lifetime distribution remains a debated issue. According to several stud- 
ies [13], species lifetime seems to be exponentially distributed. Others have found 
evidences of power law behavior with exponent close to -2 if genera, and therefore 
longer time scales, are considered ( [12], see also [11] and reference therein). Keitt 
and Stanley [14] analyzed data sets from the North American breeding bird survey 



(http://www.mbr.nbs.gov ) finding a power law distribution for species lifetime (in 
their study defined as the time between colonization and local extinction) with an 
exponent close to —3/2. In fact, the detailed analysis of Newman et al. [15] of the 
data by Raup [16] has shown how both these hypothesis consistently fit the data and, 
when a power law fit is applied, an exponent between —3/2 and —2 is estimated. 

On the theoretical side, these different, not to say contrasting, findings have been 
invoked to support different macro-ecological theories. The power law behavior with 
exponent —2 is to be expected when species dynamics can be regarded as a critical 
branching process [17] where two or more species species can originate at a random 
moment from a common ancestor and, also randomly, get extinct. An exponential 
behavior in the lifetime distribution is often referred to as Van Valen's law [18]. The 
mechanism proposed by Van Valen in support of this view is commonly known as the 
Red Queen effect: there may be no time enough for a species to gain evolutionary ad- 
vantage over competing species before the rapidly changing environment completely 
re-draws the fitness landscape. As a consequence, the extinction probability of any 
species does not depend on time and an exponential behavior for lifetimes distribu- 
tion easily follows. Several data sets support these conclusions [19-22] (see also [23] 
for further analysis of the same data). More recently, the occurrence of power law 
distributions with non trivial exponents has attracted particular attention, because 
of an ongoing debate on whether the observed patterns are caused by a self-organized 
critical dynamics [8-10] that would naturally lead to the notion of punctuated equi- 
librium [24]. In this framework an ecosystem is depicted as a system of interacting 
species whose dynamics converges spontaneously close to a critical point [12]: the 
extinction of a given species may trigger a cascade of extinction events starting from 
the species that depend upon, or directly interact with, the species just extinct and 
leading to fluctuations of any size in the number of extinction occurrences that may 
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contiguously take place. 



The aim of the present paper is to show that all the behaviors mentioned above for the 
lifetime distribution are captured by a simple model of non-interacting species. We 
conclude, therefore, that it may be problematic, if not inappropriate to discriminate 
between existing macroecological theories on the basis of existent datasets. The 
framework we adopt here is inspired by the ecological neutral theory proposed by 
Hubbell [25] and thereafter extended and analytically studied in [26-30]. This class 
of models assumes that individuals in an ecological community are fully equivalent 
and the population of a species is essentially subject to a birth and death process. 
Then, each species undergoes the same dynamics: the reproductive success of each 
individuals depends only on the species population size, and not on the particular 
species considered. Competition among species is taken into account explicitly only 
via a constraint on the total population of the community and implicitly through 
averages birth and death rates. 

From the point of view of evolutionary theory, the hypothesis of species equivalence 
may be still justified by a Red Queen effect [18], which is able to forbid the ac- 
quisition of a large evolutionary advantage (i.e. a significantly higher fitness level) 
of a species over its competitors. In the framework of population dynamics, this 
hypothesis implies that demographic stochasticity is the main driving force for the 
assembly of ecological communities, meaning that its effect is overwhelmingly strong 
compared to that of fitness differences among species, which, although present, may 
be neglected. It is worthwhile to stress that, in principle, complex ecological mech- 
anism acting on long timescales are not ruled out by these stochastic models, as far 
as they can be included in effective birth and death rates. This consideration opens 
the issue of determining whether these theories are able to assess realistic predictions 
on large time and geographical scales, such as those relevant for the fossils obser- 
vations. It is widely believed [1] that statistical physics may provide the tools to 
bridge these very different scales. In this perspective, it is encouraging that Conette 
et ai, [31] basing on the studies of the biodiversity time series compiled by Sep- 
kosky and coworkers [32], recently concluded that a random walk- like model is not 
inconsistent with the observed biodiversity time-patterns. 

The model we consider here is amenable to be analytically solved and is introduced in 
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the next section. The resulting hfetime distribution interpolates, through a scaling 
function, between the behaviors of two well known stochastic processes: exit time 
problem for the one dimensional random walk [33] and the critical Galton- Watson 
branching process [17]. Our results show that, even in a simple model in which 
interactions among species are included only in an averaged way, a variety of different 
behaviors for the distribution of extinction time is possible. In particular, depending 
on the relevant time-scales, we find an exponential, or a power law behavior. The 
latter can either occur with exponent —2, typical of branching processes [17], or, for 
shorter time-scales, with a random walk like exponent —3/2. In addition, if we assume 
that the abundance of species is distributed according to a Fisher log-series [34], in 
the Galton- Watson case, we find a power law distribution of extinction times with 
exponent —1. 

As we will discuss in the conclusion, these results stress the importance of time- 
scales and sampling effects in the analysis of lifetime distributions. This theory, 
also, can easily accommodate the contrasting empirical observations of Refs. [11- 
13, 15] by assuming that, while species lifetimes probe the exponential regime of 
the theory, genera lifetimes fall in the power law range. The fact that power laws 
arise in an "effective " single-species theory, combined with the sparseness of available 
empirical data, suggests that it may not be possible to validate (or discard) ecological 
or evolutionary mechanisms like self-organized critical dynamics [12] on the basis of 
an observed non-exponential behavior in the lifetime distributions. 



Description of the model 

According to the assumption of neutrality [25,26], the dynamics of our model is 
uniquely specified by the effective birth and death rates fe*-"^ and rf^"^ that depend 
exclusively on the population size n. 

We refer to the functions 6^*^^ and d'-"^ as "effective" because they may embody, in 
a cumulative way, a variety of ecological causes that may, in principle, influence 
the increase/decrease over time of the number of individuals in a species, or, more 
generally, in a given taxon. The framework is therefore ample enough to describe 
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a population dynamics that is not simply dominated by demographic stochasticity, 
but also, for example, by immigration, emigration or niches assembly. We can safely 
assume that h^^^ jn and S-^^n^ the birth and death rates per individuals, can be 
expanded in a power series in 1/n around their asymptotic values hi and di [28]: 

fcW/n = 61 + hQ/n + 6_i/n2 + . . . 
rf(")/n = di + do/n + d_i/r? + .... (1) 

The non-zero coefficient in this Taylor series can be generally related to various kind 
of ecological effects giving advantages (or disadvantages) to a less abundant species 
with respect to a more abundant one. In Hubbell's theory [25, 26] the terms 60 and 
df) maybe interpreted as the result of an immigration/emigration mechanism which 
couples the community to a meta-community living on a larger geographical scale. 
The mechanisms described by higher power in 1/n in Eq.(^ are relevant only for 
small population sizes and they are unable, reasonably, to affect properties observed 
on large spatial scales and long timescales. In the following, therefore, we will study 
the dynamics described only by the first two terms in the expansion of Eq. (P): 

feW = 60 + 

rfW = do + din (2) 

for all n > 1. Despite the simple form of the birth and death rates, and the simplicity 
of the assumptions, this class of models is able to provide very good fits of species 
abundance relation [26-28] which can be related to the probability -P„ of having 
species with population n. This probability, P„, evolves with time according to a 
birth and death master equation: 

We impose hi < di, ensuring that the average number of individuals is finite and 
there is no "demographic explosion". The ratio a = hi/di fixes, in-fact the average 
population per species [27,28]. 

In order to study the lifetime distribution, we consider an absorbing barrier at n = 0, 
imposing h^^^ = d^^^ = 0. The initial condition is that the new species at time t = 
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has just one individual: 

Pn{0) = Sn,l (4) 

Making these assumptions, Po{t) represents the probabihty of being aheady extinct 
at time t and the hfetime probabihty distribution function (or exit time distribution), 
p{t), is just the time derivative of PQ{t): 

Pit) = jPoit). (5) 

We win first examine the two hmit cases 61 = rfi = and = = 0^ and then move 
to the general case. 



Results 

When fei = di = the number of individuals belonging to a given species undergoes 
a random walk in n space where 60 {do) is the probability per unit time to jump one 
step to the right (left). A species lifetime would therefore correspond to the time it 
takes to the random walk to reach n = 0, i.e. to exit the positive axis. The problem 
of exit time distribution for a random walk process has been widely studied in the 
literature (see, for example, [33]). In particular, it is well known that in the critical 
case feo = do the lifetime follow a distribution of the form p{t) ~ t"^/^. Indeed it is 
easy to verify that the solution of Eq. (3), in the present case, is: 

Po(t) = 1 - exp(-2t)(/o(2t) + h{2t)) 

Pk{t)= exp(-2t)(4_i(2t)-4+i(2t)) k>Q, (6) 

where Ik{z) = ^ exp{z cos{9)) cos{k9)d9 are modified Bessel functions of integer 
order and the unit of time has been chosen such that bo = dQ = 1. Since for large z, 
Io{z) a / ^/z from Eqs. (0) and © it follows that p{t) ~ asymptotically. 

Let us now analyze the case bo = do = 0. This limiting case is interesting from an 
ecological point of view because the do and 60 terms happen to be small when one 
looks on a very large scale (like on continental scale). The dynamics is equivalent 
to a Galton- Watson process in continuous time [17]: the asymptotic behavior of the 
lifetime distribution is a classic result of the theory of critical branching processes [35] . 
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Also in this case, the birth and death equation can be analytically solved: defining 
the characteristic function G{x,t) = Yl^=o ■^ri{t)x"' the birth and death equation can 
be transformed in a first-order p.d.e. for the function G, which can be integrated 
with the characteristics method. In the following, without loss of generality, we set 
di = 1 and the initial condition in Eq. translates in G{x, 0) = x. As shown in 
details in the Supplementary Material the exact solution is: 



2 

a 



a 



This distribution has an exponential-like shape when [di — bi) (or 1 — a) is not too 
small. On the other hand, when bi approaches di, the distribution has a power law 
behavior with exponent —2 and a characteristic timescale t* = The distribution 
p{t) can be casted in a more appealing form by using the language of critical phe- 
nomena in statistical mechanics. For large t and t/t* fixed, it follows, from Eq. ([7|). 
that: 

Pit) = ^/(^) (8) 

where f{x) = [x/{l — e~^)]^e~^. Thus plotting t^p{t) versus t/t* one get, in the 
scaling region, a universal curve where all the model details are absorbed in the 
characteristic time scale, t*. When dealing with observational data, an estimate of 
t* can be obtained by the ratio of two consecutive moments, {t'') {k> 1), of lifetime 
p.d.f. . 

It is also interesting to investigate the role of the initial condition on the lifetime p.d.f.. 
Taking into account an effective speciation rate, one can show, for the particular case 
at hand, that the resulting stationary distribution [26] is the celebrated Fisher log 
series [34]: 

a" 

Pn = M— (9) 

n 

where n > and A/" is a normalization constant. Using the result above it is therefore 
possible to calculate the expected extinction time of a species that is chosen at random 
in the ecosystem. Setting as initial conditions the characteristic function associated 
to the distribution (0), G{x, 0) = log(l — xa)/ log(l — a), one finds: 



G(0,t) = — T^^g 

log(l - a) 



^l_a)e(i-")* 
e(i-a)t _ Q, 



(10) 
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In this case, again p{t) = dG{0, t)/dt e when t ^ t* whereas p{t) ~ t ^ when 
t <^ t*, which means that the critical exponent for the hfetime p.d.f is now —1. 

We now discuss, quahtatively first, the solution in the general case when all the 
coefficients are different from zero and Bq ~ do, bi ~ di. Heuristically, long-living 
species have typically a large number of individuals. For such species the bo and do 
terms can be reasonably neglected. Thus, one expects a crossover from the to the 
t^"^ behavior at a certain characteristic time and finally an exponential decay beyond 
another chi 
shown in F 
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Figure 1: numerical sample of the lifetime probability distribution function with 
parameters di = 1, bi = 1 — 5 ■ 10~^, bo = do = 10. Notice the crossover between the 
two power-laws (shown in the picture, notice the log-log scale) and the beginning of 
the exponential regime. 

In the Methods section, we provide the analytical solution of the general case, proving 
rigorously both the asymptotic critical behaviors and the scaling of the solution with 
the ratio bo/bi. In the following, instead, we discuss the main results and their 
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consequences. In terms of the Laplace transform of Po(t), Po{s) = Poit)e the 
exact solution for the critical case ho = do = r and 61 = c/i = 1, is given by 

sPo{s) - 1 = - ^<^dy,-sy^i_ly = a. log iV(.,r) ^^^^ 

where we have defined N{s,r) = ^e~*^(l — ^y. For small s the function N{s,r) 
diverges as —clogs, where c depends only on r; this implies that Po{s) behaves as 
Po{s) ~ 7 + clogs. The Tauberian theorem ensures in this case that Po{t) behaves 
like 1 — ct~^ for large t, implying that the lifetime distribution has a power-law 
tail. 

In order to derive the crossover to the behavior, we need to focus on time-scales 
t <C 1/r for r ^ 1. This is related to the limit s — with rs fixed in the solution, 
for which one obtains 

iV(s,r)=/ — e-^("+^) = 2iro(2v^) (12) 
Jo ^ 

where Ko is a modified Bessel function. Using this result and Eq. ()11|) one gets that 
the lifetime distribution obeys the following scaling form: 



Pit) = r'f (^) (13) 

where f{x) const when x — > 00, leading to the scaling at large t, and f{x) ~ \/x 
when X — > 0, corresponding to the random walk scaling at intermediate t. The 
validity of this scaling law is numerically confirmed (see Fig|2]). 



Discussion 



As sketched in the Introduction, the fact that species lifetimes are usually exponen- 
tially distributed is often referred to as Van Valen's law [18]: under the assumption 
that the fitness level is correlated in some way to the extinction probability. Van 
Valen states that an observed exponential lifetime distribution is the fingerprint of 
an acting Red Queen mechanism. Later, more detailed analysis (and datasets) [15] 
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Figure 2: Plot of curves with different values oi Bq = do = r, shown in the legend. 
The other parameters are di = 1, bi = 1 — 10"^. Collapse of the curves according to 
the scaling law (fT^ is shown in the inset. 

brought to the observation of power-law behaviors in genera lifetimes, while species 
exponential lifetime distribution have been generally confirmed. This difference is, to 
a certain degree, counterintuitive as one would expect to see a deviation from criti- 
cality as a finite size effect when looking at long time scales. A possible explanation 
proposed in Refs. [11, 12] is that at longer timescales, like those relevant for genera 
extinction, collective events like mass extinctions play a more important role. The 
interdependence of generic taxa in an ecosystem generates stronger correlations in 
their probability to survive, and these correlations, in turn, may originate a power-law 
behavior in the lifetime distribution [11,12]. 

We have shown that also in a simple model, in which every species undergoes an 
effective independent dynamics, a critical behavior for the lifetimes may occur, with 
an exponent which is compatible with the observed value. This critical behavior is 
generated only by demographic stochasticity, which is known to be a very important 
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factor in causing species extinction [36]. Interestingly enough, the hypothesis under- 
lying this model are not so different to that bringing Van Valen to the explanation 
of the exponential species lifetime: our results clearly indicate that the presence of a 
Red Queen effect, i.e. the fitness equivalence of all species, do not ensure an expo- 
nential lifetime distribution, as far as one takes into account the population sizes in 
an explicit way. In some sense, in these models the population size acts as a simple 
"memory" of the evolutionary history of the species. 

It is worthwhile to connect our approach with a model proposed by Raup [37] as a 
null model for the survivorship curves of Phanerozoic genera (the lifetime distribution 
can be thought as the derivative of the survivorship curve). This model, fitting rather 
well the fossils data, assumes that species constituting the genera have a constant 
speciation and extinction rate. Obviously, the resulting lifetimes distribution is the 
same that we recover as a limiting case of our model in Eq.©: the only difference is 
that, in Raup's case, the branching-like dynamics is applied at the level of species (not 
at the level of individuals). This implies that our model is well compatible with the 
data from the fossils record, with the advantage of being grounded on more realistic 
(and testable) hypothesis than the assumption of constant species immigration and 
speciation rates. 

In our framework, it is also possible to explain why the critical behavior in the 
lifetimes is generally observed when studying higher taxonomic levels. Let us assume 
that we can neglect the terms bo and do, as far as we are interested in the tail of the 
lifetime distribution. By taking the mean value of the distribution in Eq. Q, the 
typical population size can be expressed as: 



The r.h.s of Eq. ()14j) diverges when a 1^: thus, a choice of the parameters closer 
to criticality implies a larger population size. Since genera lump the individuals of 
may species, the effective value of a for a genera should be closer to 1 than it is for 
species. Therefore, it may not be possible to observe the power-law in the species 
lifetime due to the experimental error bars and the presence of the exponential cutoff 
occurring at t ~ (1 — a)~^ according to eq.© which, depending on its value, might 
mask both scaling regimes, i.e. and or only the latter. 




a 



1 



(14) 
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Finally, we demonstrated that, while 'local' birth and death terms, i.e. terms that are 
negligible in the large population size limit, are known to modify the mean species 
extinction time [38,39], they are unable to affect the long timescale behavior of 
the lifetime distribution: the critical behavior of the distribution, in this class of 
models, is uniquely determined by the Galton- Watson part of the dynamics. Given 
the robustness of this 'criticality' with respect to modification of the dynamics on 
small scale, we suggest the hypothesis that the observed power law could be simply 
a consequence of the branching-like structure of single population dynamics, rather 
than an effect of the interactions among different species. 



Methods 

In these notes we outline the main steps of the derivation of the results. Let us 
firstly focus on the limit = cIq = when the process reduces to a Galton - Watson 
branching process. Without loss of generality, we can set di = 1 and 6i = a in what 
follows. Introducing the characteristic function 

oo 

Gix,t) = Y,Pn{t)x\ (15) 

n=0 

the birth and death equation can be transformed in a first-order p.d.e. for G{x,t) 

dtG{x, t) = (ax^ + 1 - (a + l)x)d^G{x, t). (16) 

This equation can be integrated using, for example, the characteristic method (see 
ref. 1). Taking as initial condition G(x, 0) = x, which corresponds to Eq. 4 in the 
main text, the complete solution is 

Gix t) = nr) 

^'^ a(l-x)-(l-ax)ea-")*' ^ ' 

from which we obtain 

1 _ e(i-")* 

P(0,t) = G(0,t) = ^-^(^, (18) 

and, taking the time derivative of this, we derive Eq. 7 of the main text. It is also 
easy to see that in the scaling limit, i.e. for t* = 1/(1 — a) ^ 1 and t/t* fixed, p{t) 
can be cast in the scaling form 8. 
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In order to deal with the general case, we make a Laplace transform with respect to 
time of the generating function and define 



/>00 fOO _ _ 

G{x,s)= dte-''G{x,t) = dte"^* VP„(t)a:" 
Jo Jo 



(19) 



Then the equation of the dynamics becomes 



[ax^ + 1 — (a + l)x] dxG{x, s) + 
+ [boX + ^-bo~do-s][G{x,s)-go{s)]=sgo{s)-x, (20) 

where we defined go{s) = G{0, s), which is the Laplace transform of Poit), the function 
we wish to compute. Defining F(x, s) = G{x, s) — go{s) and using the fact that go{s) 
does not depend on x, we obtain the following equation for F{x, s): 



d^F{x, s) + p{x, s)F{x, s) = q{x, s), (21) 



where 



p{x,s) 



do bo — doa s 
X 1 — ax {1 — ax){l — x) 



( N sgo{s)-x 

= — -. (22) 

(1 — ax){l — X) 

Eq. 7 should be solved with the boundary conditions 

F(l,s) = --go{s) (23) 

s 

F(0,s) = 0. (24) 

Due to the presence of singularities at x = and x = 1, some care must be taken 
when imposing these conditions on the general solution of Eq. 7. Our strategy is 
that of solving Eq. 7 with a modified initial condition (Eq. 9) at a; = 1 — e 

F{l-e,s) = --go{s). (25) 

s 

Then we will impose condition 10 on the resulting expression, which leaves us with 
an equation for go{s). Finally, we shall restore the boundary condition 9 by taking 
the limit e ^ 0. Such an e- "regularization" procedure allows us to circumvent the 
problem of dealing with the singularities at a; = 1 of Eq. 7. Notice that, as long as a = 
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bi < 1 = di, one has \imt_,oo Po{t) = 1, i.e. the probabihty of being asjTiiptotically 
extinct approaches 1. 



The generic form of the solution of Eq. 7 with boundary condition 11 is 

^ - .9o(s) - / 



F{x, s) = e 



dx'q{x',s)e^^^ (26) 



The resulting expression is rather complex and it will be considered later on. We 
shall first specialize to the particular case Bq = = r and a = 1 discussed in the 
main text, which describes the crossover between the two power law regimes, and 
then the sub-critical case a < 1. 

For bo = do = r and a = 1, the coefficients take the simpler form 

p{x,s) = '-- jYz^ 
g(x,.) = fM_£. (27) 

Up to the leading order in e, the solution is 

F{x, s) = ^ ^ ' ■''^ ^'-'^ . (28) 

x^e 1-^ 

Since the denominator diverges when x — > 0, in order to have F{0, s) = 0, we have 
to impose that the numerator should be equal to zero. After taking the limit e — * 0, 
this yields an equation for go{s), which reads 

'^.■!S0zlf,-^^O. (29) 

Finally, upon making the substitution Y~i — V rearranging terms, we arrive at 
our main result, Eq. 11 of the main text with N{s,r) given by 



f°° dv 1 
N{s,r) = / ^e^'y{l--Y. (30) 

Ji y y 

For r fixed and s ^ 1, the integral in N{s,r) is dominated by the region y ~ 1/s 
and hence N{s, r) ~ — log s; the application of the Tauberian theorem (see ref. 2) 
finally demonstrate the asymptotic behavior of the lifetimes. In order to derive 
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Eq. 12 of the paper, in the hmit s ^ 1 with rs fixed, we make the change of variables 
X = \/^y in Eq. 16 , exponentiate the term (1 — l/yY in the integral and make a 
power expansion 

N{s,r) = r ^e-v^H^-v^^+-), (31) 



which, neglecting corrections of order \fsjr leads to Eq. 12 of the main text. When 
rs ^ 1, i.e. for t <^ r ^ 1, we can use the asymptotic expansion for the modified 
Bessel function, Kq (see Eq.l2 of the main text) or, more directly, we can estimate 
the integral with the saddle point method: the maximum of the argument of the 
exponential occurs at x* = 1 and, expanding it to second order around x* = 1, we 
find 

/■oo 

iV(s,r) ^ e-2v^ / dxe-v^(^-^)' ^ e-2^(rs)-i (32) 



Hence 

1 



s^o(s) - 1 = ^-j —. r = ^ 1 , (33) 

d, log iV(s, r) + ^ 



which means that for s — > 0, sg^i^s) — 1 ~ —\fs corresponding, according to the 
Tauberian theorem, to the random walk behavior Pi^if) ~ '^/Vt. The fact that the 
scaling variable in the derivation above is rs, implies that the crossover time should 
be proportional to r. Indeed using Eqs.ll and 12 of the main text and the inverse 
Laplace transform one derives the scaling form 

Pit) = ^/ (z) > (34) 



where the function f{x) ~ ^/x for small value of the argument (i.e. when x ^ 1) 
and approaches a constant when x becomes large. 

Finally, let us discuss the sub-critical case hi < di. Using exactly the same strategy 
as for the critical case, we find that the condition -F(0, s) = leaves us with the 
following equation: 
"1 

Citt'^« (1 - at)''o/"-rfo-^/(l-a)-l(l _ - t) = 0. (35) 







Now, we substitute y = 1 —t and solve for go{s) 

J^dy{l- yY'^ll - a{l - y)]bo/a~do~s/(l-a)-l ys/il-a) 
SQnlS) — i — ^ . (3d) 

dy (1 - yYo[i _ a{l - y)]bo/»-do^s/{i~ayi ys/ii~ayi 
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The integral on the numerator is finite when s ^ 0, whereas that on the denominator 
has a leading singularity of order (1 —a)/s. This implies that sgo{s) ^ — ^/[l + 
with A constant and t* ~ 1/(1 — a), which is exactly the Laplace transform of a 
distribution of the form 

This confirms both the asymptotic exponential decay of p{t) and the scaling of the 
cutoff time t* ~ 1/(1 - a). 
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